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Description 

[METHODS AND SYSTEMS FOR 
RESISTIVITY ANISOTROPY FORMATION 

ANALYSIS] 

Cross Reference to Related Applications 

[0001] This application claims priority pursuant to 35 U.S.C. 

§119 of U.S. Provisional Application Serial No. 60/319,479 

filed on August 19, 2002. This provisional application is 

hereby incorporated by reference in its entirety. 
Background of Invention 

[0002] Field of the Invention 

[0003] The invention relates generally to the analysis of under- 
ground earth formations, and, more particularly, to the 
determination of formation resistivity and hydrocarbon 
saturation. 

[0004] Background Art 

[0005] The localization and quantification of low resistivity pay 
(sometimes known as low contrast pay) is a primary goal 



of oil and gas exploration. Hydrocarbon saturations or 
bulk hydrocarbon volume (BVH) in subsurface isotropic 
and homogeneous rocks are usually estimated through 
the use of Archie type equations, which relate the porosity 
and water saturation in the pore space of the rock to the 
resistivity of the rock. For example, BVH in a formation 
may be determined from the total porosity and the 
total water saturation (S ) of the formation according to 

wt 

the relationship: BVH = <t> (1- S ). While the total poros- 

T wt 

ity (<fc T ) may be determined with various tools (e.g., NMR, 
neutron, density, sonic etc.), the total water saturation (S 
) often is derived from resistivity and porosity measure- 

wt 

ments, using the Archie equation or its variants (generally 
referred to as "Archie equations"). This approach of deriv- 
ing S from resistivity and porosity measurements usually 

wt 

works satisfactorily once the parameters of the equations 
are fine tuned with measurements or tests performed in 
the area. However, this approach of using Archie equa- 
tions to derive S breaks down with relation to 

wt 

anisotropic formations, and new methods are needed in 
this situation. 

[0006] Formation resistivity can be measured with conventional 
logging instruments equipped with electrodes (laterologs) 



or antennas that can transmit a current or electromagnetic 
(EM) energy into earth formations (induction/propagation 
logs). The instruments are disposed within a borehole 
traversing the formations and energy is transmitted into 
the formations to interact with the conductive media in 
the formations. With laterologs, a current (or voltage) is 
injected (or applied) into a formation using a first pair of 
electrodes and a second pair of electrodes are typically 
placed at a distance from the first pair of electrodes to 
measure the voltage drop or current flow between the 
second pair of electrodes. The measured voltage drop or 
current flow may then be used to derive the resistivity (or 
its inverse, conductivity) of the formation. 
[0007] with induction/propagation logging, EM energy is trans- 
mitted into a formation to induce eddy currents in the 
formation. The eddy currents flow in loops that lie in 
planes perpendicular to the magnetic dipole of the trans- 
mitting antenna. The magnitudes of the eddy currents de- 
pend on the conductivities of the formation. The eddy 
currents in turn induce secondary magnetic fields, the 
magnitudes of which depend on the magnitudes of the 
eddy currents. Therefore, by measuring the magnitudes of 
the secondary magnetic fields (using a receiver antenna), 



it is possible to indirectly determine the resistivity of the 
formation around the transmitter and receiver antennas. 

[0008] Resistivities of earth formations depend on the amounts 
and types of fluids included in the pores therein. Thus, 
different formations may have different resistivities due to 
different porosities, and/or different amounts or types of 
fluids included therein. When the formation is isotropic 
and homogeneous, its electric properties are constant re- 
gardless of the direction of the measurements. However, 
earth formations often comprise geological 
(sedimentation) layers that may have different petrophysi- 
cal properties (e.g., porosities, saturations, grain sizes 
etc.), and hence electrical properties. Consequently, resis- 
tivity measurements may produce different results when 
measured in different directions. This phenomenon is re- 
ferred to as formation (electrical) anisotropy. 

[0009] | n a typical situation, a borehole may be drilled through 
multiple sedimentation layers in a direction perpendicular 
to the layers, i.e., a vertical well with no formation dip. In 
such a vertical well, a resistivity measurement along a di- 
rection parallel the borehole axis is referred to as a verti- 
cal resistivity (R ) because the measurement is made in a 

V 

direction perpendicular to the sedimentation layers. In the 



vertical resistivity measurements, the current paths run 
through various sedimentation layers, which act like dif- 
ferent resistors connected in series. Therefore, the appar- 
ent resistivity of the formation (R ) is a summation of re- 

V 

sistivities contributed by individual layers. For example, as 
described below by the present invention, in a two- 
electrical-layer model comprising a relative- 
ower-resistivity layer and a relative-higher-resistivity 
layer finely interlaced together, 

R Tr =R t VF t +R hr VF hr , (1) 

where R is the measured (apparent) vertical resistivity, VF 
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and R |r are the volume fraction and resistivity of the rela- 
tive-lower-resistivity layer, and VF and R are the vol- 
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ume fraction and resistivity of the relative- 
higher-resistivity layer. As seen from Equation (1), the 
vertical resistivity measurements (R ) would be dominated 
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by more resistive layers (e.g., hydrocarbon-bearing lay- 
ers). In contrast, highly conductive thin layers (e.g., brine- 
bearing layers) may be obscured by the more resistive 
layers in such measurements. 
[0010] | n contrast, if a resistivity measurement is performed 
along a direction perpendicular to the borehole axis (or 
equivalently parallel to the beddings) in such a vertical 



well, it is referred to as horizontal resistivity (R ) because 

h 

the measurement is made along the sedimentation planes. 
Note that the "horizontal" and "vertical" used herein are 
with respect to the sedimentation layers, rather than with 
respect to the borehole axis. In horizontal resistivity mea- 
surements in a vertical well, currents flow in planes per- 
pendicular to the borehole axis, i.e., within sedimentation 
layers. Thus, each individual sedimentation layer provides 
a conductive path for the currents, i.e., the sedimentation 
layers form parallel circuits. As a result, the measured 
conductance (1/R ) is a summation of the conductivities 

h 

of the sedimentation layers within the region of investiga- 
tion. For example, in the electrical-layered model com- 
prising relative-lower-resistivity and relative- 
higher-resistivity layers, 

1 PR VF. 



where R h is the horizontal resistivity, VF and R |r represent 
the volume fraction and resistivity of the relative- 
lower-resistivity layer, and VF and R represent the vol- 
ume fraction and resistivity of the relative- 
higher-resistivity layer, respectively. As seen from Equa- 
tion (2), the measured horizontal conductivity (R ) would 

h 



be dominated by the most conductive layers, while thin 
non-conductive hydrocarbon-bearing layers may become 
"invisible" to the horizontal resistivity measurements. This 
is known as the low resistivity pay problem. 

[0011] As mentioned above, either the conductive layers or the 
resistive layers could be "missed" in the resistivity mea- 
surements depending on the direction of the measure- 
ments, in the sense that their presence has only a small 
effect on the apparent resistivity. Thus, formation resistiv- 
ity anisotropy presents a problem in formation evaluation. 
Many reservoir rocks exhibit resistivity anisotropy. Several 
mechanisms can produce the anisotropy, among which 
are very thin sand (carbonate)-shale laminations, grain 
size changes in clean sandstone, wind-distributed sands 
(aeolian formations), cementing (porosity) changes in 
sandstone and so on. See Rubin, D. M., Cross bedding, bed- 
forms, and paleocurrents, Society Of Economic Paleontolo- 
gists And Mineralogists, Concepts In Sedimentology And 
Paleontology, 1; Klein et al., The petrophysics of electrically 
anisotropic reservoirs, Transactions of the SPWLA Thirty- 
Sixth Annual Logging Symposium, Paris, France, June 
26-29, 1995, paper HH. 

[0012] over the years, most of the homogeneous or thick-layer 



oil and gas reservoirs have been discovered. As a result, 
many reservoirs comprise thin layers of pay zones. With 
technology advances such as directional and horizontal 
drilling, it is becoming economical to produce in thin 
reservoirs that traditionally would have been ignored. The 
industry has also begun to realize the importance of thinly 
laminated reservoirs that have been by-passed due to low 
apparent resistivity in vertical wells. Therefore, a need ex- 
ists for methods that can accurately predict the hydrocar- 
bon contents of thinly laminated reservoirs. 
[0013] Evaluation of thinly laminated reservoirs is not a new 
problem in formation evaluation and interpretation. See 
U.S. Pat. Nos. 3,166,709 and 5,461,562, assigned to the 
present assignee. The horizontal and vertical resistivities 
of anisotropic formations can be evaluated by wireline or 
logging-while-drilling (LWD) EM measurements in highly 
deviated wells. See Hagiwara T., A New Method to Determine 
Horizontal-Resistivity in Anisotropic Formations with Prior Knowl- 
edge of Relative Dip, Transactions of the SPWLA Thirty- 
Seventh Annual Logging Symposium, June 16-19, 1996, 
Paper Q; U.S. Patent Nos. 5,966,013, 6,092,024, 
5,886,526, 6,218,841. Recent techniques using tri-axial 
EM tools and the combination of EM tools with laterolog 



tools have made it possible to measure the vertical and 
horizontal resistivities in vertical wells. However, it re- 
mains a difficult task to relate the measured resistivities 
(R and R ) to the bulk hydrocarbon volume because the 

v h 

traditional Archie relation does not apply to anisotropic 
formations. 

[0014] Several papers have been published dealing with analyses 
of anisotropic formations. See Klein et al., The Petrophysics 
of Electronically Anisotropic Reservoirs, Transactions of the SP- 
WLA Thirty-Sixth Annual Logging Symposium, June 
26-29, 1995, Paper HH; Tabanou et al., Which Resistivity 
Should Be Used To Evaluate Thinly Bedded Reservoirs at High An- 
gle? 'Transactions of the SPWLA Fortieth Annual Logging 
Symposium, May 30- June 3, 1999, Paper E; Shray F. and 
Borbas T., Evaluation of Laminated Formations Using Nuclear 
Resonance and Resistivity Anisotropy Measurements, SPE Eastern 
Regional Meeting, Canton, Ohio 17-19 October 2001; U.S. 
Pat. No. 5,550,473. 

[0015] These papers propose methods for determining electrical 
properties of anisotropic formations. As shown in Equa- 
tions (1) and (2), four parameters (VF , R , VF , and R ) 

hr hr Ir Ir 

are determined in order to define the electrical properties 
of a two-layer model. With 



VF„ + VF k =l, (3) 

there are four unknowns and three equations. Note that R h 
and R are assumed known parameters that can be deter- 
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mined from resistivity measurements. Thus, one addi- 
tional parameter is needed to solve these equations. 
[0016] Depending on the additional parameter that is used to 
solve these equations (hence, electrical properties of the 
formation layers), conventional methods may be catego- 
rized as the 

method, the 

method, and the 

R^-R^ -BFV 

method, where R , V , and BFV are shale resistivity, shale 

sh sh 

content, and bound fluid volume, respectively. Once the 
electrical properties of the formation layers (e.g., R and 
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R ) are known, they can be used together with the total 
porosity (<t> T ) of the formation to determine the total water 
saturation (S ). The total water saturation (S ) and total 

wt wt 

porosity (4> T ) may then be used to determine the bulk hy- 
drocarbon volume (BVH) according to the relationship BVH 



= c(> (1-S ). It is assumed that a reliable technique exists 

T wt 

for the water saturation evaluation of each of the individ- 
ual homogeneous layers in the layered model. The Archie 
relation is typically assumed applicable in the individual 
layers. 

[0017] The additional parameter used in these methods is either 
obtained from other types of measurements (e.g., NMR, 
GR) or assumed to be the same as the value determined in 
a thick layer of the same composition. However, these as- 
sumed values may not accurately represent the values of 
the same types of layers in thin laminations. In addition, 
these methods are based on assumptions of formation 
geological compositions and petrophysical properties (i.e., 
grain sizes or porosities). If the formation has a different 
composition or property from that assumed, these meth- 
ods cannot provide accurate estimates of electrical prop- 
erties of the formation layers. As a result, the derived S 

wt 

and BVH may not be accurate. 
[0018] Thus a need remains for improved techniques that can 
provide accurate S and BVH estimates without prior 

wt 

knowledge of the properties of the thin layers of the for- 
mation. 
Summary of Invention 



[0019] one aspect of the invention provides a method for deter- 
mining a subsurface formation property. The method in- 
cludes determining a volume fraction of a layer in a multi- 
electrical-layer model for an anisotropic region of the for- 
mation, wherein the multi-electrical-layer model includes 
a relative-lower-resistivity layer and a relative- 
higher-resistivity layer and the determining is based on a 
resistivity measurement; determining a resistivity for the 
relative-lower-resistivity layer and a resistivity for the rel- 
ative-higher-resistivity layer based on the volume fraction 
and bulk resistivity measurements of the anisotropic re- 
gion, wherein the bulk resistivity measurements include a 
horizontal resistivity measurement or a vertical resistivity 
measurement; and determining the formation property 
based on the volume fraction, the resistivity of the rela- 
tive-lower-resistivity layer, the resistivity of the relative- 
higher-resistivity layer, a porosity measurement of the 
anisotropic region, and the bulk resistivity measurements. 

[0020] Another aspect of the invention provides a system for de- 
termining a subsurface formation property. The system 
includes a computer system adapted to process a pro- 
gram including instructions for: determining a volume 
fraction of a layer in a multi-electrical-layer model for an 



anisotropic region of the formation, wherein the multi- 
electrical-layer model comprises a relative- 
lower-resistivity layer and a relative-higher-resistivity 
layer and the determining is based on a resistivity mea- 
surement; determining a resistivity for the relative- 
lower-resistivity layer and a resistivity for the relative- 
higher-resistivity layer based on the volume fraction and 
bulk resistivity measurements of the anisotropic region, 
wherein the bulk resistivity measurements include a hori- 
zontal resistivity measurement or a vertical resistivity 
measurement; and determining the formation property 
based on the volume fraction, the resistivity of the rela- 
tive-lower-resistivity layer, the resistivity of the relative- 
higher-resistivity layer, a porosity measurement of the 
anisotropic region, and the bulk resistivity measurements. 
[0021] other aspects and advantages of the invention will be ap- 
parent from the following description and the appended 
claims. 

Brief Description of Drawings 

[0022] FIG. 1 shows geological models which can cause forma- 
tion (electrical) anisotropy and their equivalent model 
(multi-electrical-layer model) in accord with the invention. 

[0023] FIG. 2 shows a high-resolution resistivity measurement 



used to determine a volume fraction of a relative- 
lower-resistivity layer in accord with the invention. 
[0024] FIG. 3 shows a flow chart of a process in accord with the 
invention. 

[0025] FIG. 4 shows high-resolution and low-resolution resistiv- 
ity measurements of a formation comprising isotropic re- 
gions and anisotropic regions in accord with the inven- 
tion. 

[0026] FIG. 5 shows a conventional computer system that may be 

used to implement a system according to the invention. 
Detailed Description 

[0027] Embodiments of the invention relate to general methods 
and systems for determining a formation property such as 
a total water saturation (S ) from bulk resistivity mea- 

wt 

surements (R and R ) in anisotropic formations. The total 

v h 

water saturation (S ) together with the total porosity (cb ) 

wt T 

may then be used to determine the bulk hydrocarbon vol- 
ume (BVH). The disclosed techniques use formation mod- 
els having layers with different electrical properties. These 
formation models are referred to as "electrical models," as 
opposed to the geological models of known techniques. 
The term "electrical" is used to indicate that the formation 
model corresponds to the electrical characteristics of the 



formation as opposed to a conventional geological model. 
Conceptually, the electrical model is based on the inher- 
ent electrical properties of the formation regardless of the 
geological properties. Although the electrical properties of 
a given layer are dependent on the composition of the 
formation, one need not know the makeup of the layer 
because the resistivity measurements sense the electrical 
properties of the formation rather than the chemical or 
physical compositions of the formation. 
[0028] while embodiments of the invention use electrical models 
that are independent of actual formation compositions, 
conventional methods often rely on models that are spe- 
cific to formations composed of particular compositions. 
For example, the 

method assumes a two-geological-layer model compris- 
ing a sand layer and a shale layer. The sand layer com- 
prises coarser grains and includes hydrocarbons in the 
pores, hence a higher resistivity. On the other hand, the 
shale layer comprises finer grains and often includes wa- 
ter in the pores, hence a lower resistivity. 
[0029] | n t he 



method, the additional parameter used to determine the 
resistivities of the two layers is the thin shale layer resis- 
tivity (R ), which is assumed to be identical to the re- 
shale 

sistivity (R ) determined in the nearby thick shale layer. 

sh 

The assumed R value (R in Equations (1) and (2)) to- 
sh Ir 

aether with R,R are then used to compute R (i.e., R 

v h hr sand 

in this case) and VF (i.e., VF in this case). This model 

Ir shale 

further assumes that the porosity of the sand layer is the 
same as that of the shale layer (hence, 



) in order to determine S from the computed R and R 

wt sand 

. The estimated S is then used to determine bulk hy- 

shale wt 

drocarbon volume (BVH) according to the relationship: 

bvh = m-^j- 

BVH = (1-S ). 
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[0030] This method is accurate if the anisotropic formation in- 
deed includes only shale and sand layers, i.e., this method 
depends on whether the studied formation is equivalent to 
the formation described in the sand-shale two-ge- 
ological-layer model. Further, the R estimate obtained 

3 1 shale 

from nearby thick shale may not represent the true value 
of R in the thin layers, and the assumption that the 

shale 

porosities of the sand layer and the shale layer are identi- 



cal may not be valid. Furthermore, this method impliedly 
assumes R > R , an assumption that may not be 

sand shale 

correct in some formations. 
[0031] Similarly, the 

method assumes that the volume fraction of the thin shale 
layers (VF ) is identical to the shale content (V ), i.e., 

shale sh 

VF = V . Based on this assumed VF (i.e., VF in 

shale sh shale Ir 

Equations (1) and (2)) and the R and R measurements, R 
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(i.e., R in this case) and R (i.e., R in this case) are 

sand Ir shale 

then computed by solving Equations (1) (3). Next, this 
method assumes the porosity of the sand layer is identical 
to that of the shale layer 

in order to calculate S from the determined R and R . 

wt hr Ir 

The calculated S is then used to derive BVH. 

wt 

[0032] Like the 

method, the 

method is typically accurate for anisotropic regions in- 
cluding only shale and sand layers. Furthermore, the shale 



content (V ) may not correctly represent the volume frac- 

sh 

tion of shale layers (VF ) if the shale layers include 

shale 

coarse grain non-clay-minerals such as quartz, and the 
sand layers include clay. This method also impliedly as- 
sumes R > R ,an assumption that may not be cor- 

sand shale 

rect in some formations. 

[0033] The 

R^-R^ -BFV 

method uses the bound fluid volume (BFV) determined 
with other sensors (e.g., NMR) and the total porosity (d> T ) 
to derive the volume fraction of the thin fine grain layers 

(VF ), i.e., 

fg 

VF fg = BFV0 r . 

With VF (i.e., VF in Equations (1) and (2)) known, R 

fg Ir hr 

(i.e., the resistivity of coarse grain layers R in this case) 
and R (i.e., the resistivity of the fine grain layers R in 

Ir fg 

this case) can then be determined. This method further 
assumes that the porosity of the coarse grain layer is 
identical to the porosity of the fine grain layer 

Based on this assumption, S is then calculated according 
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to an equation similar to the Archie equation. The calcu- 



lated S is then used to derive BVH. 
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[0034] T he 

R^-R^ -BFV 

method is accurate for anisotropic formations including 
only shale and sand, or for formations in which one layer 
is a coarse grain clean reservoir with lower water satura- 
tion and the other layer is fine grain clean formation con- 
taining no hydrocarbons and both layers have the same 
porosity. Again, this method depends on whether the 
studied formation is equivalent to the formation described 
in the sand-shale two-geological-layer model or in the 
coarse-fine grain two-geological-layer model. Further- 
more, the assumption that 

VF fg =BFV/js5 r 

may not correctly represent the actual value of the volume 
fraction of fine grain (low-resistivity) layers, depending on 
whether the fine grain layers include free fluid or the 
coarse grain (high-resistivity) layers include bound water. 
This method also assumes 

= K = & " 

and impliedly assumes "R > R ", assumptions which 

eg fg 

may not be correct for some anisotropic formations. 



[0035] | n contrast to these conventional methods, embodiments 
of the invention use "electrical models" of formation lay- 
ers. These electrical models do not rely on prior knowl- 
edge of the formation layer compositions. Therefore, em- 
bodiments of the invention may be generally applied in all 
formation types, regardless of the underlying mechanism 
(saturation changes, porosity changes or both) for the 
anisotropy. 

[0036] a multi-electrical-layer model of the invention entails dif- 
ferent types of "electrical" layers, including one having a 
relative higher resistivity and the other having a relative 
lower resistivity. Note that the multi-electrical-layer 
model refers to differing types of "electrical" layers, which 
may include more than two types of "geological" layers. 
The difference in resistivity between these electrical layers 
may arise from different porosities, different grain sizes, 
different formation fluids and/or different water satura- 
tions of the various types of geological layers. 

[0037] An anisotropic formation according to a multi-electri- 
cal-layer model embodiment of the invention has bulk re- 
sistivity measurements (R and R ) as defined in Equations 
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(1) and (2). Together with Equation (3), there are three 
equations and four unknowns. Therefore, one additional 



parameter is determined before the resistivities (R and R 

hr 

) of the multiple layers are solved. 
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[0038] | n accordance with embodiments of the invention, the ad- 
ditional parameter may be a volume fraction of the rela- 
tive-lower-resistivity layers (VF^) or the relative- 
higher-resistivity layers (VF^). While embodiments of the 
invention may use either VF or VF , for clarity the fol- 
lowing description refers to VF^ only. In one embodiment 
of the invention, VF^ may be determined from a high- 
resolution resistivity measurement. The high-resolution 
resistivity measurement may be acquired with any con- 
ventional logging tool known in the art. Some of these 
tools do not provide accurate (absolute) resistivity mea- 
surements. Instead, they provide relative resistivities of 
the various layers in the formation. Embodiments of the 
invention may use high-resolution resistivity measure- 
ments that include either absolute or relative resistivity 
values. Note that in accordance with the invention, the VF 

Ir 

is determined from actual electrical measurements rather 
than using an assumed value or a value determined from 
non-electrical measurements as in prior methods. 
[0039] | n accordance with an embodiment of the invention, VF 

Ir 

may be determined from the high-resolution resistivity 



measurement by finding the bed boundaries of the vari- 
ous layers and then assuming the thicknesses of the low- 
resistivity layers within an index interval (sliding window). 
The sliding window can be any size, but preferably larger 
than the vertical resolution provided by the resistivity 
measurement. A default size of the sliding window may be 
equal to the vertical resolution of a low-resolution resis- 
tivity tool which provides bulk resistivity measurements 
containing Rv and Rh signals. 
[0040] Finding the bed boundaries from the high-resolution re- 
sistivity measurement may be accomplished by taking a 
derivative of a curve representing the high-resolution re- 
sistivity measurement as a function of the measurement 
depth (or borehole axial depth). The bed boundaries 
would appear as inflection points on the derivative curve. 
Once these inflection points are identified, the alternate 
layers defined by these points represent layers having the 
same resistivity (e.g., relative lower resistivity). Thus, the 
thicknesses of these layers may be summed to provide a 
volume fraction of the relative-lower-resistivity layers, VF 
, for example. 

Ir 

[0041] FIG. 1 shows various geological models which can cause 
formation anisotropy and their equivalent multi-elec- 



trical-layer model according to embodiments of the in- 
vention. Since a source of the apparent formation 
anisotropy is that the thickness of any one of a group of 
(at least two) adjacent geological layers with different re- 
sistivity values is lower than the vertical resolution obtain- 
able from the resistivity measurement within this group of 
layers, the multi-electrical-layer model can represent all 
the multi-geological layer models which may cause for- 
mation anisotropy. The difference between resistivity val- 
ues among geological layers can be caused by different 
porosities or different saturations. In reality both porosity 
and saturation can be different in each of the layers. 

[0042] FIG. 2 shows a high-resolution resistivity measurement 
that is used to determine bed boundaries in accordance 
with an embodiment of the invention. As shown, a high- 
resolution resistivity log (400) includes a high-resolution 
resistivity curve (401) that represents the high-resolution 
resistivity measurement plotted against the y-axis (404) 
that indicates the borehole axial depth at which the resis- 
tivity was measured. The x-axis (402) in this plot indi- 
cates relative resistivities. 

[0043] Also shown on the high-resolution resistivity curve (401) 
are inflection points 411-418 as determined by taking 



derivatives of the high-resolution resistivity curve (401). 
As noted above, these inflection points 411-418 repre- 
sent the bed boundaries between layers of different resis- 
tivities. As shown in FIG. 2, these inflection points (bed 
boundaries) define three relative-lower-resistivity layers 
A, B, C interleaved with relative-higher-resistivity layers 1, 
2, 3, 4. According to this example, VF |r may be calculated 
by the sum of the thicknesses of layers A, B, and C divided 
by the sum of the thicknesses of layer A, B, C, 1, 2, 3 and 
4. 

[0044] | n an alternative embodiment, the bed boundaries may be 
approximated as points (not shown) defining the full 
width half maximum (FWHM) of peaks on the high- 
resolution resistivity curve (401). The thicknesses of the 
alternating layers of relative-lower-resistivity defined by 
these points are then summed to give the VF |r in a manner 
similar to that described for the inflections points 
411-418. 

[0045] | n y e t another embodiment, a cutoff curve (406) may be 
calculated based on a selected criterion, such as a curve 
linking the inflection points 411-418 or the FWHM points. 
Alternatively, the cutoff curve (406) may be a low- 
resolution resistivity curve acquired with a low-resolution 



tool and having its magnitudes adjusted to be comparable 
to that of the high-resolution resistivity curve (401). With 
the cutoff curve (406) defining layers of different resistivi- 
ties, VF may be calculated as: 

lr 



w ir ■ (4) 



[0046] | n one embodiment of the invention, 

corresponds to the sum of the thicknesses of all the layers 
where the high-resolution resistivity curve (401) is lower 
than the cutoff curve (406) within an index interval 
(sliding window). Accordingly, TH corresponds to the 

sw 

thickness of all the layers within an index interval (sliding 
window), i.e., the thickness of the sliding window. 
[0047] Once the VF is determined, the resistivities of the two 

lr 

layers, R and R , may be determined by solving Equa- 
tions (1) (3). In one embodiment, Equations (1) (3) may be 
rearranged to give Equations (5) and (6) to facilitate this 
calculation. 
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[0048] The calculated values of R and R are then combined 
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with 



(the total porosity) to compute BVH and S using the 

wt 

equations similar to the Archie-type equations, for exam- 
ple, Equations (7) and (8): 



BVH = $ T -[(I - R?) + R? J|l ] (7) 



^-i ^ . (8) 

[0049] Note that in these computations, only the total porosity ^ 
is used. There is no need to separate the porosities of the 
relative-higher-resistivity and relative-lower-resistivity 
layers, <b and <b . That is, embodiments of the invention 
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do not rely on porosities of the individual layers and, 
therefore, no error is introduced by improper assump- 
tions. The total porosity 4> t may be obtained from con- 
ventional logging operations as known in the art. 
[0050] Equations (7) and (8) assume that the layers contain one 
type of water, or the free water resistivity is more or less 
equal to the clay bound water resistivity. In this case, the 
Archie equation can be used (according to the Archie 



equation, 



, where, ^ and 



is the water porosity and electrical porosity of the forma- 
tion, R is the water resistivity, and R is the true resistiv- 

w t 

ity of the formation). If the free water resistivity is not 
more or less equal to the clay bound water resistivity, the 
dual-water model may be more appropriate. According to 
the dual-water model, 



R.. 



R, 



(?) 



where 

$w and ikUe 

is the water porosity and electrical porosity of the forma- 
tion, R is the effective water resistivity, and R is the 

we t 

true resistivity of the formation. Equation (9) is analogous 
to the Archie equation. The effective water conductivity 
(1/R or c ) may be defined as: 

we we 



1 



R^ L - 



S . 



c_ + 



(10) 



where S is the total water saturation, s is the clay 

wt wcl 



water saturation, c is the free water conductivity, C is 

w clw 

the clay bound water conductivity, and VF is the volume 

clw 

fraction of clay bound water, c may be determined us- 

we 

ing nuclear magnetic resonance (NMR) data as known in 
the art. For example, C = (1 VF ) C + VF C , 

we clw w clw clw 

where VF = BFV / BFV, BFV is clay bound water vol- 

clw cl cl 7 

ume, BFV\s clay and capillary bound water volume. C 

w 

may be determined from a water reservoir, and C may 

clw 

be computed from pure shale. 
[0051] Once the effective water conductivity (C or 1/R ) is 

we we 

determined, the total water saturation (S ) may be ob- 

wt 

tained as: 



s . = 1 - — 



■Hi 



(11) 

[0052] Finally, the bulk hydrocarbon volume (BVH) may be com- 
puted from the total porosity (<)> T )and the total water satu- 
ration (S ): BVH = 4 (1-S ). 

wt T wt 

[0053] FIG. 3 shows a flow chart illustrating a process of the in- 
vention. Initially, a volume fraction of a layer in a multi- 
electrical-layer model is determined for an anisotropic re- 
gion of the formation, wherein the multi-electrical-layer 
model includes a relative-lower-resistivity layer and a rel- 



ative-higher-resistivity layer and the determining is based 
on a resistivity measurement (Step 300). A resistivity for 
the relative-lower-resistivity layer and a resistivity for the 
relative-higher-resistivity layer are also determined based 
on the volume fraction and bulk resistivity measurements 
of the anisotropic region, wherein the bulk resistivity 
measurements include a horizontal resistivity measure- 
ment or a vertical resistivity measurement (Step 305). The 
formation property is then determined based on the vol- 
ume fraction, the resistivity of the relative- 
lower-resistivity layer, the resistivity of the relative- 
higher-resistivity layer, a porosity measurement of the 
anisotropic region, and the bulk resistivity measurements 
(Step 310). 

[0054] The bulk resistivity measurements include vertical resis- 
tivity (R ) measurements and horizontal resistivity (R ) 

v h 

measurements acquired using any tool capable of provid- 
ing the R and R resistivities as known in the art, includ- 

v h 

ing conventional EM induction or propagation resistivity 
tools. 

[0055] a high-resolution resistivity measurement is preferred in 
Step 300. Those skilled in the art will appreciate that em- 
bodiments of the invention may use any conventional 



high-resolution resistivity measurement tools, e.g., tools 
sold under the trade names of FMI™, RAB/GVR™, SHDT™, 
OBMI™, OBDT™ by Schlumberger Technology Corporation 
(Houston, TX). For example, if wireline tools are used with 
a water based mud, then FMI™ or SHDT™ may be used to 
obtain a high-resolution resistivity measurement. If wire- 
line tools are used with an oil based mud, then OBMI™ or 
OBDT™ may be used to obtain a high-resolution resistivity 
measurement. If LWD tools are used in a water-based 
mud, then RAB™ may be used to obtain a high-resolution 
resistivity measurement. If LWD tools are used in an oil 
based mud, then a micro-induction type of high- 
resolution resistivity sensor such as the OBDT™ sensor 
may be used to obtain a high-resolution resistivity mea- 
surement. 

[0056] N 0 te that whether the high-resolution is acquired before, 
at the same time, or after the acquisition of the bulk re- 
sistivity measurements is immaterial to embodiments of 
the invention. The high-resolution resistivity measure- 
ment is usually not calibrated, and thus the resistivity val- 
ues are only meaningful in a relative sense. Even after cal- 
ibration, these high-resolution resistivity curves often do 
not represent true formation resistivities because the 



high-resolution resistivity measurements are typically ac- 
quired in regions close to the borehole and are likely con- 
taminated by the drilling fluids. Because the high- 
resolution resistivity curve is only used to calculate a vol- 
ume fraction of the relative-lower-resistivity layers (or the 
relative-higher-resistivity layers) according to embodi- 
ments of the invention, it is not critical that the high- 
resolution resistivity measurement reflect true formation 
resistivities or relative resistivities, i.e., accurate resistivity 
values are not required to implement the invention. 
[0057] Embodiments use the resistivity measurement to derive a 
volume fraction of the relative-lower-resistivity layers, VF 
, (or the relative-higher-resistivity layers, VF ). As noted 
above, VF^ may be derived using one of several methods. 
For example, the bed boundaries between layers of differ- 
ent resistivities may be identified as inflection points on a 
derivative curve of the high-resolution resistivity curve. 
The thicknesses of the relative-lower-resistivity layers are 
then summed and divided by the total bed thickness to 
give the VF . Alternatively, a cutoff curve may be deter- 
mined and used to derive the VF |r , for example, according 
to Equation (4). 

[0058] The described measurements and the VF determined are 

Ir 



used to calculate the electrical properties (e.g., resistivi- 
ties, R and R ) of both the relative-lower-resistivity and 

Ir hr 

the relative-higher-resistivity layers. The R ^ and R may 
be determined according to Equations (5) and (6) de- 
scribed above. R and R are determined with actual elec- 

lr hr 

trical logging data (high-resolution electrical measure- 
ments and R and R ), and no assumption or estimate of 

v h 

an additional parameter is required in deriving these val- 
ues. Thus, the obtained R and R values are more reli- 

lr hr 

able and no error is introduced by improper assumptions 
or improper geological models. 
[0059] The R and R hr together with the total porosity (<j>^ are 

used to derive the total water saturation, S , and the bulk 

wt 

hydrocarbon volume (BVH). These calculations may be ac- 
complished as shown in Equations (7) and (8) or Equation 
(11), depending on the water model. Again, note that the 
S and BVH are determined from the obtained data and 

wt 

the derived R and R ; no assumed value is required. 

Ir hr' M 

Thus the techniques of the invention provide more reliable 
determination of these parameters. 
[0060] Those skilled in the art will appreciate that the techniques 
of the invention may be applied to formations including 
isotropic and anisotropic regions. For example, referring 



to FIG. 4, formation (500) includes four regions (A, B, C, 
and D). Each region includes one or more layers. Further, 
regions A, C, and D are isotropic, while region B includes 
thin layers that have different resistivites. Accordingly, 
high-resolution resistivity measurements (502) of regions 
A, C, and D do not exhibit significant variance from the 
corresponding low resistivity measurements (504). How- 
ever, the high-resolution resistivity measurements (502) 
of region B shows significant variations from the corre- 
sponding low-resistivity measurements (504). The exis- 
tence of variance between the high-resolution resistivity 
measurements (502) and the low-resolution resistivity 
measurements (504) indicates that the region is 
anisotropic. 

[0061] Techniques of the invention involve determining VF by 
locating bed boundaries from the high-resolution resis- 
tivity measurements. These methods are not aimed at 
identifying boundaries in regions A, C, and D shown in 
FIG. 4, but will find bed boundaries in Region B. Thus the 
techniques of the invention may be applied to any forma- 
tion in general, with formations including anisotropic re- 
gions being resolved into layers having different electrical 
properties and isotropic regions not being affected. 



[0062] The above description is based on electric models having 
two types of formation layers having different resistivities. 
Note that the two "type" of electrical layers may include 
more than two types of geological layers. Depending on 
the resolution of the high resolution resistivity measure- 
ment, one may dissect the anisotropic region into smaller 
regions and use the multi-electrical-layer models dis- 
closed herein to solve the electrical properties of each 
smaller region. 

[0063] Some embodiments of the invention relate to systems for 
performing the techniques described above. A system of 
the invention may be a stand-alone unit for performing 
methods of the invention or may be incorporated in a 
drilling tool (i.e., an onboard processor unit). A system in 
accordance with the invention typically includes a proces- 
sor and a memory. In some embodiments, a system may 
be implemented on a general-purpose computer having a 
processor, a memory, and may optionally include other 
hardware. For example, as shown in FIG. 5, a typical com- 
puter (50) includes a processor (52), a random access 
memory (54), and a storage device (e.g., permanent mem- 
ory or hard disk) (56). The computer (50) may also include 
input means, such as a keyboard (58) and a mouse (60), 



and output means, such as a monitor (62). Note that the 
general-purpose computer is only for illustration and em- 
bodiments of the invention may take other forms (e.g., in- 
tegrated in a logging tool, not shown). Note also this pro- 
cess can be carried out during or after the data acquisi- 
tion. In case of while-drilling applications, the process can 
be either performed on surface with data sent uphole in 
real-time via telemetry and/or downhole on acquired data 
and then sending the processed results to surface via 
telemetry. Note also that this process can be implemented 
as part of a general formation evaluation program or as an 
independent evaluation program by itself. 
[0064] | n a system in accordance with the invention, the memory 
stores a program readable by the processor. The program 
may include instructions for performing the above de- 
scribed methods, e.g., determining a volume fraction for a 
layer in an anisotropic formation based on a high- 
resolution resistivity measurement; determining formation 
electrical properties based on the volume fraction and 
bulk resistivity measurements that include horizontal re- 
sistivity and vertical resistivity measurements; calculating 
formation total water saturation based on the electrical 
properties of the formation; and determining a bulk hy- 



drocarbon volume of the formation. 

[0065] a system in accordance with the invention provides new 
and improved techniques to evaluate the properties of 
anisotropic reservoirs in real time. The programming may 
be accomplished through the use of one or more program 
storage devices readable by the computer processor and 
encoding one or more programs of instructions exe- 
cutable by the computer for performing the operations 
described herein. The program storage device may take 
the form of, for example, one or more floppy disks; a CD- 
ROM or other optical disk; a magnetic tape; a read-only 
memory chip (ROM); and other forms of the kind well 
known in the art. The program of instructions may be in 
"object code," i.e., in binary form that is executable more- 
or-less directly by the computer; in "source code" that re- 
quires compilation or interpretation before execution; or 
in some intermediate form such as partially compiled 
code. The precise forms of the program storage device 
and of the encoding of instructions are immaterial here. 

[0066] The present invention provides several advantages. Meth- 
ods of the invention provide a general approach to cor- 
rectly translate R and R into S for anisotropic forma- 

v h wt 

tions. The methods are based on multi-resistivity-layer 



electrical models which are compatible with multi- 
geological-layer formations. Methods of the invention fo- 
cus on the real cause of resistivity anisotropy: the differ- 
ent resistivity values, rather than the different formation 
compositions. Accordingly, the disclosed technique sim- 
plifies a wide range of geological models into a single cal- 
culation model (i.e., an "electrical" model) that may simply 
include a relative-higher-resistivity layer (R ) and a rela- 

hr 

tive-lower-resistivity layer(R ). These methods do not rely 

Ir 

on the assumption of "4> = <t> ." The methods in- 

K ^lr ^hr T 

elude high-resolution resistivity measurements to derive a 
volume fraction of a layer; no assumed value is used to 
solve the electrical properties of the formation layers. Fur- 
ther, the techniques may be used to obtain continuous R 

Ir 

and R curves without problems because R is always 
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greater than R |r according to the electrical models. 
[0067] while the invention has been described with respect to a 
limited number of embodiments, those skilled in the art, 
having benefit of this disclosure, will appreciate that other 
embodiments can be devised which do not depart from 
the scope of the invention as disclosed herein. The 
present invention may be implemented using conventional 
EM logging instruments and in general oilfield operations 



such as while-drilling, wireline applications, coiled tubing 
operations, while-tripping, monitoring applications, etc. It 
will also be appreciated that the present invention may be 
implemented in conjunction with other models to deter- 
mine subsurface properties as known in the art. 



